Microrheological Characterisation of Anisotropic Materials 
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We describe the measurement of anisotropic viscoelastic moduli in complex soft materials, such 
as biopolymer gels, via video particle tracking microrheology of colloid tracer particles. The use of a 
correlation tensor to find the axes of maximum anisotropy without prior knowledge, and hence the 
mechanical director, is described. The moduli of an aligned DNA gel are reported, as an application 
of the technique; this may have implications for high DNA concentrations in vivo. We also discuss 
the errors in microrheological measurement, and describe the use of frequency space filtering to 
improve displacement resolution, and hence probe these typically high modulus materials. 

PACS numbers: 83.80.Lz, 83.10.Mj, 87.15. La, 87.14.Gg, 83.80.Xz, 83.60.Bc 



Many important biological molecules are capable of 
forming liquid crystalline phases due to their inherent 
stiffness and aspect ratio. In this paper we demon- 
strate how a development of existing microrheological 
techniques offers a means to characterize local anisotropic 
properties in such media. 

We believe that this technique is unique in its abil- 
ity to determine the mechanical properties of typical 
anisotropic materials, where the director may vary over 
a mesoscopic lengthscale. Our approach also utilises sev- 
eral improvements to existing microrheology measure- 
ments, which are likely to be particularly important for 
such high modulus, anisotropic materials for which there 
are particular challenges. Although we have chosen to 
work with DNA, many materials, both naturally occur- 
ing and synthetic, form anisotropic phases. Examples of 
such materials of biological importance include aqueous 
biopolymer gels, including DNA, but also cellulose, xan- 
than and F-actin, amongst others. Anisotropy may also 
occur in intracellular materials, to which standard 2-D 
microrheology has previously been applied [lj, d ■ 

Particle tracking microrheology is a technique that has 
attracted much interest in experimental soft condensed 
matter as a method for measuring the viscoelastic proper- 
ties of soft materials on the micron scale. The mechanical 
response of such a material may be probed by monitor- 
ing the time evolution of the displacement of embedded 
microscopic tracer particles subject to thermal motion. 

The idea of using embedded particles to probe ma- 
terial properties was first described in a seminal paper 
by Freundlich and Seifriz Q, studying gelatin gels using 
magnetic beads in an external field. However these ex- 
periments suffered from the use of aspherical particles, 
limited force control and inaccuracy in the detected par- 
ticle position. The technique was rediscovered by Mason 
and Weitz |5j who proposed an extension to the familiar 
Stokes-Einstein diffusion equation to complex, frequency 
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where the real viscoelastic shear modulus G of a material 
at a temperature T may be calculated as a function of 
Laplace frequency, s, for beads of radius a, with a dis- 
placement r, in d dimensions, given the mean squared dis- 
placement (MSD), (r 2 (r)) = (|r(i + r) - r(*)H. In ad- 
dition to improved experimental techniques, this presen- 
tation of a generalised Stokes-Einstein relation (GSER) 
has encouraged renewed interest, and resulted in many 
advances, greatly expanding the power of microrheology 
as a technique (for a review see e.g. @). 

From the MSD, under the same assumptions as the 
generalised Stokes-Einstein relation (GSER) (EqnQ we 
can calculate the creep compliance by multiplication 
with a constant, 
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or extract the frequency dependent complex viscoelastic 
shear modulus via a Laplace or Fourier transform 0, 
Eo| . If the material is a gel, i.e. a material composed of 
a polymer network in a fluid medium, then the network 
spacing (mesh size) may also be determined via particle 
tracking [TT| . 

Conventional microrheology has so far used the 2- 
dimensional MSD in Eq. QJ, but we choose to decom- 
pose the displacement into 2 orthogonal 1-dimensional 
displacements. Although an individual random walk is 
spatially anisotropic |l2(, the probability density func- 
tion averaged over a number of walks will be circularly 
symmetric in an isotropic material. In a mechanically 
anisotropic material, the ensemble average of the ran- 
dom walks will also be elliptically symmetric (Fig. ^) 
with major and minor axes along the mechanical direc- 
tor. The GSER thus allows us to calculate the dissipa- 
tive and elastic properties of the gel along specific axes, 
subject to the same conditions as the Stokes derivation: 
no-slip boundary conditions on a sphere in an incom- 
pressible continuum fluid, with no inertia. In addition 
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FIG. 1: (a) Measured probability density plot of probe po- 
sitions in the aligned DNA gel described later (left) and a 
similar unsheared DNA sample (right), as an example of an 
isotropic material. The plots are calculated by collapsing 
probe tracks (longer than 50 steps) over all available time 
to a common origin, and counting the occurrences in discrete 
bins. The angle of the surrounding ellipse is set at 35° from 
the horizontal, which is the value calculated from the eigen- 
vector analysis for the shortest lag time later. A good corre- 
spondence is seen, (b) Displacement distribution in x (left) 
and y (right) for r = l/25s. Visually the diffusion is Gaus- 
sian, and this may be verified by calculating the excess kur- 
tosis of the distributions (0.8 in the s-direction and 0.7 in the 
y- direct ion ) . Thus the degree of heterogeneity is relatively 
small, and may be neglected. 



the implicit assumption is made that the Stokes drag for 
viscous fluids may be generalised to viscoelastic materi- 
als at all frequencies, and that the probe particles have a 
negligible effect on the material 0, ^1 . Although in gen- 
eral viscosity is an 81 element 4th rank tensor, many of 
these elements are either zero or repeated in typical cases. 
For example a nematic liquid crystal has 5 independent 
intrinsic viscosities describing diffusion. Our case is is 
analogous to the dra g on a s phe re moving though a ne- 
matic liquid crystal [Til ITU llr| and so, since we have 
a rotational symmetry axis parallel to the director, the 
Brownian motion is governed by two independent diffu- 
sion coefficients where D±/\\ = k^T/ (6irr]x/\\)a, as the 
derivation of the Stokes-Einstcin relation is separable in 
spatial dimension. Generalising to complex viscosity, we 
thus expect the GSER to hold in anisotropic materials, 
with two independent complex viscosities r]^_ and 77jj . The 
ratio, rf\_/rfn is therefore a measure of the anisotropy of 
the system. Similarly to the isotropic case, it should be 
noted that there is no rigorous theoretical backing for 
the extension to complex moduli, only phenomeological 
justification. 

The probe particles themselves may be tracked via a 
number of methods. We choose real-space multiparticle 
'single-point' video tracking due to both the large amount 
of information available for analysis, and the relative ease 
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FIG. 2: Mean squared displacement of 0.497/^m polystyrene 
latex colloidal spheres dried on a glass slide, calculated from 
all independent lag times. The deinterlaced video fields ap- 
pear displaced from each other by ~ 0.15 pixels (most clearly 
seen at short lag times), varying by ~ 0.06 pixels with time 
(for clarity even lag times, in frames, are denoted by empty 
symbols) . This represents a significant error for high modulus 
materials. 

of setup. Video is acquired with a Zeiss-Axioplan opti- 
cal microscope with a lOOx oil immersion lens (N.A. = 
1.30) and conventional commercial CCD video camera 
(Sony SCC-DC138P PAL), directly into a computer us- 
ing a Hauppauge Impact VCB capture card. The analysis 
is performed using custom-written MATLAB scripts. To 
detect the small Brownian displacements of colloidal par- 
ticles in high modulus materials we need as high a spatial 
resolution as possible. This may be measured by deter- 
mining the apparent displacement of fixed spheres, fol- 
lowing [TtJ ■ We use a 2-dimensional bandpass filter to se- 
lect appropriate spatial frequencies in Fourier space, and 
thus filter out noise and background variation. Regions 
above an intensity threshold are selected, then checks for 
eccentricity and shape solidity are performed to deter- 
mine the quality of the probe images. Finally the cen- 
troid of the region is calculated, weighted by the original 
pixel intensity. By this method we have achieved a ID 
displacement resolution of up to 3nm, with 4nm in typi- 
cal operation. 

Since the video data in conventional cameras is 
recorded in an interlaced manner, in which the two fields 
of alternate lines are recorded either 1/50 or 1/60 of a 
second apart in each frame, one typically de-interlaces 
the fields and analyses each field separately. However, 
by tracking fixed colloidal spheres at high resolution, we 
show in Fig. that the two fields are displaced from 
each other in a time varying fashion. The source of this 
displacement is unclear but it means that de-interlacing 
may be inadvisable when tracking small displacements, 
and is not carried out in the current work. 

High modulus gels often have some in-built shear stress 
when they are made, which relaxes over a very long time 
period. In addition one may wish to investigate the mi- 
crorheological properties of sheared, flowing, or even liv- 
ing samples. However, probe particles in such a system 
will be subject to a net flow field, in addition to Brow- 
nian motion. Provided that the flow field is uniform, 
(Ar(r)) = for a random walk, so the net drift velocity 



can be calculated from the ensemble average and then 
subtracted. If the flow field varies across the field of view 
this approach is clearly invalid, and if the velocityfield is 
a function of depth the situation is complicated [18j , but 
in practice we find this simple approach to be adequate in 
this work, where we measure the drift due to relaxation 
as ~ 3/im/min. 

We write the displacement in tensor form to remove 
coordinate system dependence from the arbitrary x (hor- 
izontal) y (vertical) Cartesian axes with the origin in 
the top-left in the video. The eigenvector relation (Eqn. 
|3J) can then be written, and the resulting eigenvectors 
give the basis in which displacements along the two or- 
thogonal axes are least correlated, and hence maximally 
anisotropic. We use the eigenvectors calculated from the 
shortest lag time, r, for which there are the most events 
and hence the smallest error, to find the new basis and 
hence recalculate displacements in these coordinates. For 
an anisotropic material one would expect the new axes to 
lie parallel and perpendicular to the average optical di- 
rector. By solving the equation for the specific example 
shown in Fig. |3 and determining the eigenvectors (su- 
perimposed on Fig. |3J, we can indeed see the correlation 
of the optical alignment with the directions of anisotropy. 
The directionality and strength of alignment in the mate- 
rial can thus be determined from particle tracking alone, 
with no prior knowledge required. We will discuss the 
nature of the sample below. 
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The MSD must then be computed. This is a con- 
cise representation of the rheological data, related to the 
creep compliance by a constant (Eqn. |2J) . We need to de- 
termine the errors in microrheological data, with a partic- 
ular view to judging whether anisotropy is significant, so 
the following expression is derived, making the assump- 
tion that a 'typical' value for x is \J (Ax 2 (t)) . Assuming 
Gaussian statistics, we have that a x i = 2a x x and so, 
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where N is the number of events contributing to 
(Aie 2 (t)). The expression is of the appropriate form as 
<j x = y/ (Ax 2 (t)) in this case. In addition, the noise at 
the resolution limit results in a 'static error' 01 which 
is exhibited as a constant offset added to (Ar 2 (r)), and 
thus can be subtracted. 

A number of methods have been proposed to arrive at 
a complex shear modulus 0, M, OH D3 • We use a simple 
discrete Laplace transform, 
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FIG. 3: The sheared llmg/ml DNA sample at 298K under 
horizontal-vertical crossed polarizers. Alignment is relatively 
weak. Light areas show alignment, with the DNA chain back- 
bones lying approximately 45° counterclockwise from the ver- 
tical (checked via insertion of a quarter waveplate at 45° to the 
polarizers). The calculated eigenvectors are superimposed in 
black, and show the mechanical alignment corresponds well 
with the optical birefringence, as expected. The calculated 
angle is 35° counterclockwise. 



where N is the number of samples spread over the interval 
t = to t = tjv, followed by 1st order polynomial spline 
fitting and analytic continuation into the complex plane, 
substituting iu> for s in the locally fitted form |5|, |2(| • No 
assumptions are made about the form of the response 
which is well approximated by a linear spline fit over a 
5 point window, and truncation errors are found to be 
small for a discrete Laplace transform given a typical 
range of data in time (l/25s to 60s). 

To demonstrate the utility of microrheology as a tool 
for probing anisotropic materials, we apply the method 
described above to shear-aligned deoxyribonucleic acid 
(DNA). DNA in water, with some ammonium acetate as 
a counter ion, forms well-documented liquid crystalline 
phases 0j at high concentrations. We look at genomic 
DNA (Fluka, Deoxyribonucleic acid sodium salt from 
herring testes. Product: 31162) at a moderate concen- 
tration (llmg/ml) at 298K in 0.25M ammonium acetate, 
contained in a cell made of a metal washer (0.80mm 
deep) sandwiched by a glass microscope slide and cov- 
erslip, holding 12/il of sample. The sample is left to relax 
for 1 hour, then shear stress is applied by moving the 
coverslip, to induce alignment, which may be monitored 
by optical birefringence (Fig. |3J). We use 0.3/xm diame- 
ter probes (polystyrene latex, Agar Scientific) which are 
much larger than the expected mesh size of 15nm (cal- 
culated assuming uniform volume filling), with our field 
of view set well away from the cell surfaces. The video is 
kept interlaced, as discussed previously, and 120 seconds 
is analyzed. 

Although the sample exhibited only weak optical align- 
ment, we find significant mechanical differences between 
the axes parallel and perpendicular to the optical direc- 
tor. Figure 0] plots the MSD, showing greater motion 
perpendicular to the optical director than parallel. This 
seems reasonable in a network of stretched chains, with 
the optical director lying along the chain backbone. The 
viscoelastic moduli are also calculated (Fig. [5J. The 
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FIG. 4: 1-dimensional MSDs along calculated eigenvectors for 
sheared llmg/ml DNA at 298K. Crosses (+) correspond to 
displacements along the axis approximately parallel to the op- 
tical director, and circles (o) to displacements perpendicular 
to this. The static resolution error is taken as 4nm, and sub- 
tracted. The MSD is proportional to the creep compliance 0. 
Inset: MSD for the same DNA at 4.9 mg/ml, where there was 
no visible birefringence. There is now no visible anisotropy 
in the MSD. Note that the shape of the response is similar 
to the anisotropic case - this suggests that the GSER should 
hold in anisotropic materials. 



material is clearly stiffer to shear deformation in the di- 
rection along which the chains are, on average, already 
stretched. This is as we would expect for a polymer net- 
work which has been stretched uniaxially, where a pro- 
portion of chains are close to fully extended between en- 
tanglements/ crosslinks in this direction. It is interesting 
to note that there appears to be a cross-over in the elas- 
tic modulus, where the stretched direction has a lower 
elastic modulus than the perpendicular direction at fre- 
quencies below 3Hz, the explanation for which is unclear. 
The plateau in G' is typical of a gel in the elastic region of 
its frequency response. Some studies [22, UM have been 
done on DNA at lower concentrations, which indicate 
precursors to this elastic regime. 

We have presented a method to look at the viscoelas- 
tic properties of anisotropic materials on typical (meso- 
scopic) lengthscales. The technique both permits the 
identification of the principal axes of the system, and 
the associated mechanical anisotropy to be characterized 
as a function of frequency. Some rheological data for me- 
chanically aligned DNA is presented. These observations 
may have implications for high DNA concentrations in 
vivo. Several improvements to microrheology measure- 
ments have been described, with a particular focus on 



high modulus anisotropic materials. Many intracellular 
materials appear anisotropic, and the approach presented 
allows the elucidation of the mechanical properties of the 
materials, potentially in vivo |24j . Higher modulus ma- 
terials should be able to be explored using fluorescent 
probes with diameters in the 10-100nm range. It should 
be possible to extend the technique to '2-point' [2^| cross- 
correlation techniques to ensure that probe effects on the 
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FIG. 5: 1-dimensional elastic modulus, G'(ui) and loss modu- 
lus, G"(u>), along calculated eigenvectors for sheared llmg/ml 
DNA at 298K. Filled circles correspond to G' along the 
axis approximately parallel to the optical director, and filled 
squares to G' perpendicular to this. Open symbols refer to 
G"(ui) in these directions. The static resolution error is taken 
as 4nm, and subtracted from the MSD before calculation. In- 
set: Moduli of the same DNA at 4.9 mg/ml (approximately 
the mean concentration of DNA in a human cell |22j|1. The 
sample is isotropic in its mechanical response, and the shape 
is comparable to the anisotropic case, suggesting the validity 
of the GSER in anisotropic materials. 



network |2|f (27| can be neglected |13|, |28| . We anticipate 
that the technique described will also be of use in deter- 
mining the anisotropic pore sizes of materials which form 
an anisotropic network. 
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